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SUMMARY 


The report examines several buckling analysis procedures for stiffened panels, presents 
accurate results for seven stiffened panels, and illustrates buckling modes with plots of buck- 
ling mode shapes. All panels are rectangular and have stiffeners in one direction - down the 
length of the panel. 

The buckling analyses used in PASCO are summarized with emphasis placed on the 
shear buckling analyses. PASCO buckling analyses include the basic VIPASA analysis, which 
is essentially exact for longitudinal and transverse loads, and a smeared stiffener solution 
(equivalent orthotropic plate solution) that was added in an attempt to alleviate a shortcoming 
in the VIPASA analysis - underestimation of the shear buckling load for modes having a buck- 
ling half-wavelength equal to the panel length. Such buckling modes are sometimes referred 
to as overall modes or general instability modes. Buckling results are then presented for 
seven stiffened panels loaded by combinations of longitudinal compression and shear. The 
buckling results were obtained with the PASCO, EAL, and STAGS computer programs. The 
EAL and STAGS solutions were obtained with a fine finite element mesh and are very accurate. 
These finite element solutions together with the PASCO results for pure longitudinal compres- 
sion provide benchmark calculations to evaluate other analysis procedures. For each 
example, several figures illustrate buckling mode shapes for pure compression and pure shear 
loadings. It was concluded that the smeared stiffener solution should be used only with 
caution. 


INTRODUCTION 

Although buckling analysis procedures that are both fast and accurate have been devel- 
oped for stiffened panels subjected to longitudinal (N x ) and transverse (N y ) loadings (for 
example, VIPASA, refs. 1 to 3, and BUCLASP2, ref. 4), no such procedure has been developed 
for analyzing stiffened panels subjected to loadings involving shear (Njjy). (See fig. 1.) 
VIPASA nearly meets the dual objectives of speed and accuracy; however, when the loading 
involves shear, VIPASA underestimates the buckling load for the overall mode - that is, the 
mode for which the buckling half-wavelength in the direction of the stiffeners is equal to the 
panel length. VIPASA is generally accurate for loadings involving shear when the buckling 
half-wavelength in the direction of the stiffeners is less than one-third the panel length. 

Shear buckling analysis procedures in current use include the following modeling 
approaches: stiffeners modeled as linked plates with infinite panel length (VIPASA, ref. 1); 

hinges along plate element connections for local buckling and smeared stiffnesses for overall 
buckling (for example, ref. 5); approximations in which stiffeners are modeled as discrete 
lines of bending (El) and twisting (GJ) stiffnesses on the panel skin; and general purpose finite 
element approaches (for example, EAL, refs. 6 and 7, and STAGS, refs. 8 and 9). All these 




approaches have shortcomings. The shortcoming of the approach used in VIP ASA is mentioned 
in the previous paragraph and is discussed in this report. When stiffnesses are smeared, 
local deformations that contribute to the overall buckling mode are lost. Some local deforma- 
tions are also lost when the stiffeners are modeled as El and GJ stiffeners. Finite ele- 
ment approaches can provide high accuracy by using detailed modeling and fine meshes; how- 
ever, to obtain accurate results, the computation costs may be high. 

Because of the shortcomings in VIP ASA, the buckling analysis in PASCO (refs. 10 to 13), 
an alternate solution approach for predicting overall shear buckling was explored and incorpo- 
rated in PASCO. That approach is based on smeared orthotropic stiffnesses. 

This report presents buckling results obtained with the computer programs PASCO 
(which includes both VIP ASA and a smeared orthotropic solution), EAL, and STAGS for seven 
stiffened panels. For each panel, results are presented for several combinations of inplane 
shear and longitudinal (stiffener direction) compression ranging from pure shear to pure 
longitudinal compression. The results serve three purposes. They help evaluate the shear 
buckling analyses in PASCO, they provide accurate benchmark calculations to evaluate other 
analysis procedures, and they help provide a better understanding of the buckling mechanism 
for stiffened panels through the numerous detailed plots of buckling mode shapes. 

SYMBOLS 

Values are given in both SI and U.S. Customary Units. The calculations were made in 
U.S. Customary Units. 

b plate width 

D plate bending stiffness 
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orthotropic stiffnesses defined in equation (2) 


i 




D 11’ D 22’ D 33 
E Young's modulus 

®1’®'2 Young's modulus of composite material in fiber direction and transverse to fiber 
direction, respectively 

E 12 shear modulus of composite material in coordinate system defined by fiber 

direction 


El 

e 

GJ 

L 

My 

N x 

N x cr 

N xy 

N xy C r 

N y 

p 


bending stiffness of beam 

amplitude of bow-type imperfection at panel midlength 
twisting stiffness of beam 
panel length 

bending moment about line parallel to Y-axis (see fig. 2) 
bending moment about line parallel to X-axis 

applied longitudinal compressive loading per unit length (see fig. 2) 

value of N x that causes buckling 

applied shear loading per unit length (see fig. 2) 

value of that causes buckling 

applied transverse loading per unit length (see fig. 2) 

lateral pressure 


u,v,w buckling displacements 

W panel width 

X,Y,Z coordinate axes in longitudinal, transverse, and lateral directions, respectively 

x,y,z coordinates 
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strain in x-direction 


e x 

0 fiber orientation angle (see fig. 9) 

X buckling half-wavelength 

\i Poisson's ratio 

of composite material in coordinate system defined by fiber 
E 1 

= M 2Ei 

Anti indicates that buckling displacement w is antisymmetric with respect to a point 

at the center of the panel 

FACTOR eigenvalue, the product of FACTOR and applied load is buckling load 
S.S. simple support boundary conditions 

Sym indicates that buckling displacement w is symmetric with respect to a point at 

the center of the panel 

Computer programs: 

BUCLASP2 BUC kling of LAminated Stiffened Panels 
EAL Engineering Analysis Language 

PASCO Panel Analysis and Sizing COde 

STAGS Structural Analysis of General Shells 

VIP ASA Vibration and Instability of Plate Assemblies including Shear and Anisotropy 

BUCKLING ANALYSIS IN PASCO FOR LOADINGS INVOLVING SHEAR 

PASCO is a computer program for analyzing and sizing uniaxially stiffened composite 
panels subject to the loadings shown in figure 2. PASCO is described in references 10 to 13. 




Poisson's ratios 
direction, j 


Abbreviations: 
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X 
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Figure 2.- Stiffened panel with initial bow, applied 
loading, and coordinate system. 

As pointed out in the introduction, an important limitation of PASCO is that VIP ASA 
(the buckling analysis in PASCO) underestimates the buckling load when the loading involves 
shear and the buckle mode is a general or overall skewed mode having a longitudinal buckle 
length equal to the length of the panel. That limitation and an approximate analysis technique 
intended to overcome it are discussed in this section. 

VEP ASA Buckling Analysis 

VIPASA, the buckling analysis program incorporated in PASCO, treats an arbitrary 
assemblage of plate elements with each plate element i loaded by (N x )j, (Ny)j, and (N^j. 
The buckling analysis connects the individual plate elements and maintains continuity of the 
buckle pattern across the intersection of neighboring plate elements. The buckling displace- 
ment w assumed in VIPASA for each plate element is of the form 

w = f x (y) cos ™ - f 2 (y) sin H (1) 

where X is the buckling half-wavelength. Similar expressions are assumed for the inplane 
displacements u and v. Because the buckling displacements are assumed to have a speci- 
fied form in the x-direction, the VIPASA solution is essentially a one-dimensional solution. 
(The finite element solutions discussed subsequently are two-dimensional.) The functions 
fj(y) and f 2 (y) satisfy the differential equation of equilibrium and allow various boundary 
conditions to be prescribed on the lateral edges of the panel. Boundary conditions cannot be 
prescribed on the ends of the panel. However, certain useful boundary conditions are implic- 
itly satisfied at the ends of the panel. VIPASA is, therefore, still effective for analyzing a 
broad spectrum of structural analysis problems. Boundary conditions at the ends of the panel 
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and the effect of these boundary conditions on the predicted buckling load are discussed in the 
following paragraphs. 

For orthotropic plate elements with no shear loading, the solution given by equation (1) 
involves a series of node lines that are straight, perpendicular to the longitudinal panel axis, 
and spaced X apart, as shown in figure 3. Along each of these node lines, the buckling dis- 
placements satisfy simple support conditions. Therefore, for values of X given by 
X = L, L/2, L/3, L/m, where m is an integer and L is the panel length, the nodal pat- 
tern shown in figure 3 satisfies simple support boundary conditions at the ends of a finite, 
rectangular, stiffened panel. For this case, the VIPASA solution is exact. 


Node lines 

Ny 

t ♦ ♦ t t ♦ t 


N 


x 




♦ 


♦ ♦ 


♦ 

K 


♦ I 



♦ 


Figure 3.- Node lines for orthotropic plate elements 
with no shear loading. 


For anisotropic plate elements and/ or plate elements with a shear loading, the solution 
given by equation (1) involves node lines that are skewed and not straight, but the node lines 
are still spaced X apart, as shown in figure 4. (Because anisotropy generally has negligible 
effect on buckling loads for long- wavelength buckling modes and because it is these long- 
wavelength modes that are troublesome, reference to anisotropy is dropped in the following 
discussion.) Since node lines cannot coincide with the ends of the rectangular panel, the 
VIPASA solution for loadings involving shear is accurate only when many buckles form along 
the panel length, in which case boundary conditions at the ends are not important. An example 
in which X = L/4 is shown in figure 5. 

As X approaches L, the VIPASA buckling analysis for a panel loaded by may 

underestimate the buckling load substantially. One explanation is as follows. As seen in fig- 
ure 5, the skewed nodal lines given by VIPASA in the case of shear do not coincide with the end 
edges. Forcing node lines to coincide with the end edges produces buckling loads that are 
always higher than those determined by VIPASA. When only one buckle forms along the panel 
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Node lines 


"y 

t * * * * 



v 


Figure 4.- Node lines for anisotropic plate elements and/or 
plate elements with shear loading. 



Figure 5-- Buckling of panel under shear loading. 
Mode shown is X = L/M . 


length, the effect of boundary conditions at the ends is much more important than when multiple 
buckles form along the length. For long-wavelength buckling modes, the buckling load for a 
panel satisfying the end boundary conditions can be more than twice the buckling load for a 
panel not satisfying the end boundary conditions. Such cases are illustrated in the examples 
discussed subsequently. 
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In summary, for stiffened panels composed of orthotropic plate elements with no shear 
loading, the VIP ASA solution is exact in the sense that it is the exact solution of the plate 
equations satisfying the Kir chhoff- Love hypothesis. However, for stiffened panels having a 
shear loading, the VIP ASA solution can be very conservative for the case X = L. 

Because VTPASA is overly conservative in the case of long-wavelength buckling if a 
shear load is present, another easily adaptable analysis procedure based on smeared ortho- 
tropic stiffnesses has been examined for the case X = L. That attempt to get an improved 
solution is discussed in the following section. 

Smeared Stiffener Solution 

The objective of the analysis is to solve the shear buckling problem for the finite panel 
illustrated in figure 6. For buckling half-wavelength X equal to panel length L, the mathe- 
matical model solved by VIP ASA and the resulting node lines are similar to those illustrated 
in figure 7. The panel in figure 7 is infinitely long in the x-direction, but X is finite. 



Figure 6.- Finite stiffened panel of length L and 
width W, simply supported on all four edges and 
subjected to shear load N X y* S.S. designates 
simple support boundary conditions. 


It is assumed that a better approximation to the solution for the finite panel would be 
obtained with the infinitely wide panel shown in figure 8. That assumption is based on the 
belief that, for this example, it is more important to satisfy boundary conditions along the 
edges perpendicular to the stiffeners than along the edges parallel to the stiffeners. There 
are two reasons for this belief. First, stiffeners tend to cause node lines to be nearly parallel 
to the stiffeners. This means that boundary conditions along edges parallel to the stiffeners 
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are more nearly satisfied automatically. Second, stiffeners tend to cause a buckle pattern 
having more half-waves transverse to the stiffeners than in the stiffener direction. Such a 
buckle pattern causes boundary conditions along edges parallel to the stiffeners to be 
relatively unimportant. 

Unfortunately, the mathematical model illustrated in figure 8 cannot be analyzed with 
VIP ASA because VIP ASA requires that the panel be uniform in the direction of the infinite 
dimension. However, the mathematical model obtained by smearing the stiffnesses of the 
stiffened panel of figure 8 can be analyzed with VIPASA. That solution is referred to as the 
smeared stiffener solution. It is obtained by calculating orthotropic stiffnesses for the stiff- 
ened panel and interchanging the x and y loading and stiffnesses. Directions are changed 
because the x-direction is the infinite direction in VIPASA. The eigenvalue used is the lowest 
of the set for X = W, W/2, W/3, ..., where W is the panel width. 

The orthotropic stiffnesses involved in the smeared stiffener solution are the bending 
stiffnesses D-q and D22 and the twisting stiffness D33. These stiffnesses are defined in 
the following differential equation for lateral deflection of an orthotropic plate with lateral 
loading q: 


D 117T +4D 33 

8x* 


8x 2 8y 


4 

, T-V 8*W _ 

2 22 ^4 5 


( 2 ) 


The procedure used to calculate the orthotropic stiffnesses is described in reference 11. 

The attempt to improve on the VIPASA solution for long-wavelength shear buckling is 
more involved than the discussion presented here. However, the basic feature - smeared 
stiffener solution - of that solution approach and the conclusions regarding its suitability are 
the same as those presented here. A more complete discussion of the procedure is presented 
in references 11 and 13. As is shown in the following examples, the smeared stiffener solu- 
tion does, in some cases, provide an improved solution; however, the procedure should be 
used only with caution. 


STIFFENED PANEL EXAMPLES 

Seven stiffened panels were analyzed with PASCO and with the finite element structural 
analysis codes EAL (refs. 6 and 7) and STAGS (refs. 8 and 9). Results of these analyses are 
presented in this section. Four of the seven panels had blade stiffeners, one panel had hat 
stiffeners, one panel was a corrugated panel, and one panel had Z stiffeners. Four panels 
were modeled with a graphite- epoxy composite material; three were modeled with metal 
(aluminum). All panels were 30.0 in. (76.2 cm) square and had six equally spaced stiffeners. 
The loadings were combinations of longitudinal compression (N x ) and shear (N^). PASCO 
and EAL results are presented for each loading. STAGS results are presented only for the 
pure shear loadings and only for the first four examples - the blade- stiffened panels. The 
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VIP ASA results provide an accurate check of the EAL model for pure longitudinal compres- 
sion, and the STAGS results for the first four examples provide an independent check of the 
EAL results for pure shear. All standard VIPASA solutions (not the smeared stiffener 
solutions) include the effect of anisotropic bending stiffness terms for each plate element 
making up the panel cross section. The EAL and STAGS solutions also include anisotropic 
effects. 

A schematic drawing showing the loading and overall dimensions for the seven example 
cases is shown in figure 9. The direction shown for N X y in figure 9 is positive for EAL and 
STAGS (and for this report) and corresponds to a negative value of for PASCO. The 

positive orientation for the fiber orientation angle 6 is also shown. For all panels and load- 
ings, the stress analysis in PASCO (ref. 11) is used to distribute the applied loads over the 
panel cross section to obtain the prebuckling stress state. In particular, the N x load is dis- 
tributed by assuming uniform strain e x of the panel cross section with free transverse 
expansion of each plate element, so that Ny for each plate element is zero. PASCO, EAL, 
and STAGS buckling analyses all use the same prebuckling load distribution. A table contain- 
ing the prebuckling load distribution for a typical loading is given for each example. 



Figure 9.- Loading and dimensions for stiffened panel examples. 


The goal is for buckling boundary conditions to be simple support on all four edges. 
These bolmdary conditions are applied to the stiffeners as well as to the skin. The simple 
support boundary conditions can be defined as follows: 

At x = 0, L; Mx = w = 6N X = 6v = 0 
and at y = 0, W; M y = w = 6N y = 6u = 0 
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where 6N X , 6v, 6Ny, and 6u represent changes in N x , v, Ny, and u from the values 
in the prebuckling stress state. Buckling boundary conditions used by the analyses are as 
follows: 

PASCO 

At x = 0, L; boundary conditions cannot be prescribed. Simple support boundary 
conditions exist at the ends x = 0 and x = L only for the pure 
compression case. See previous discussion in section entitled 
"VEP ASA Buckling Analysis." 


<< 

u 

o 

y 

3 

V 

and 

&w 

9y 

are unrestrained, and 

u = w = 0 

EAL and STAGS 






At x = 0, L; 

u 

and 

dw 

Bx 

are unrestrained, and 

o 

II 

II 

& 

II 

> 

At y = 0, W; 

V 

and 

Bw 

By 

are unrestrained, and 

— -S-® 


In each example, the solution denoted "VIP ASA" is the standard VIP ASA solution based 
on a detailed model of the panel cross section. The solution denoted "smeared stiffener 
solution" is also calculated with VIP ASA, but the panel is modeled as an orthotropic plate. 

The panel cross sections are treated as collections of lines with no offsets to account 
for thicknesses. (Offsets are available in PASCO, but were not used in order to allow direct 
comparison with other methods that cannot easily account for offsets.) The first example is 
discussed in greater detail than the other examples. 

Example 1 - Composite Blade-Stiffened Panel 

Panel description .- A repeating element of the composite blade -stiffened panel is shown 
in figure 10. Element widths are also shown. The wall construction for each plate element is 
given in table I. Only half the laminate is defined for each plate element because all laminates 
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Figure 10.- Repeating element for example 1, composite 
blade-stiffened panel. 



are symmetric. Plate element numbers are indicated by the circled numbers in figure 10. 
Fiber orientation angles are measured with respect to the X-axis, which is parallel to the 
stiffener direction. 

Values of Young's moduli, shear modulus, and Poisson's ratio for the graphite- epoxy 
material used in the calculations for this example are given in table II. 

Table HI contains the prebuckling load distribution for the applied loading 
N x = = 1000 lb/in. (175.1 kN/m). Since the prebuckling stress analysis treats longitudi- 

nal compression and shear separately, the prebuckling load distribution for other combina- 
tions of N x and can be obtained by separately scaling the values of longitudinal com- 

pression and shear given in table HI. The orthotropic stiffnesses used for the smeared 
stiffener solution are = 81 614 lb-in. (9221.1 N-m), D22 = 2893.1 lb-in. (326.88 N-m), 
and D33 = 1601.7 lb-in. (180.97 N-m). 

PASCO input .- Sample PASCO input for this example is shown in figure 11. The input 
for this and the other examples is in U.S. Customary Units. The loading is N x = 1000, 

Nxy = -1000, which means that a solution is sought for the case N x = -N X y = 1000 lb/in. The 


***** EXAMPLE 1, COMPOSITE BLADE- STIFFENED PANEL ***** 
$C0NDAT 
$ 

$PANEL 

B=2 .5 , 1.352, 2.5, 

T= .0055 , .0055, .0495, .0055, .011, 

THET=45, 0, 90, 45, 0, 

KWALL( 1 , 1)= 1 1 ,- 1 , 1 , 2 , 3 , 

KWALL( 1 ,2)=4 ,-4,-4, 4, 5 , 

IWALIf= 1,2,1, 

HCARD=4,-4,2,90,0, 

2,121,4, 

4.5.1.3, -121, 

ICARD=5, 1,3, 1,-909, 0900, 

3. 2. 3. 4, 

3, 3, 4, 3, 

3.4, -909,0900, 

N0BAY- 6 , 

ICREP=6 , 

EIf=30 , 

MINLAM=30, 

NLAM=1,2, 

SHEAR= 1 . , 

IBC= 1 , 

IP=2 , 

NX-1000. , 

NXY--1000 . , 

$ 

$ MATER 

El-19. E6, E2-1.89E6, E12-.93E6, ANU1-.38, RHO-.0571, 

$ 


Figure 11.- Sample PASCO input for example 1, 
composite blade-stiffened panel. 
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negative sign on causes the direction of in PASCO to agree with that shown in 

figure 9. The unknown in the PASCO eigenvalue analysis is a scalar factor which multiplies 
the load vector. ICARD input is included in order to get detailed plots of buckling mode 
shapes. The repeating element shown in figure 10 was generated with PASCO input. 

EAL model .- The single finite element type used in the EAL model for this and the 
other examples is a four-node, quadrilateral, combined membrane and bending element. Both 
the membrane and bending stiffness matrices for the element are based on the assumed 
stress, hybrid formulation of Pian (refs. 6, 14, and 15). For this element, the geometric 
stiffness matrix that is required in the buckling analysis is based on a conventional displace- 
ment formulation that includes terms allowing inplane (u and v) as well as out-of-plane (w) 
displacements. The Pian membrane formulation allows a single element across the depth of 
a blade stiffener to represent accurately its overall inplane bending behavior. (Two elements 
were used to account for out-of-plane deformations.) The EAL designation for this element 
is E43. 

The finite element grid chosen for the EAL model is shown in figure 12. Two elements 
are used along the depth of the blade, 4 elements are used between blades, and 36 elements 
are used along the length. Thus, a total of 1296 elements and 1369 nodes are present. A 
more refined model consisting of 3 elements along the depth of the blade, 4 elements between 
blades, and 48 elements along the length was used as a check for selected loadings. In all 
cases, the differences in the results for the two models were negligible. Also, since the 
PASCO solution is essentially exact for the pure compression case, PASCO results served as 
a guide when developing finite element models for EAL and STAGS. 

It is recognized that the EAL model just described and the STAGS models described in 
the next section are more refined than models used in usual engineering calculations. How- 
ever, since accuracy was an issue in this study, it was decided that benchmark calculations 
that differ from the exact solution by no more than approximately 1 percent were needed. 

Based on convergence studies and other comparisons, it is believed that the finite element 
calculations presented in this report meet this accuracy requirement. 

As pointed out in a previous section, instead of solving the prebuckling stress problem 
with EAL and STAGS, the prebuckling load distribution is calculated with PASCO and is input 
to EAL and STAGS so that all three buckling analyses use a common prebuckling load distri- 
bution. Buckling boundary conditions used by EAL and STAGS were given on page 12. Addi- 
tional details of the EAL input and analysis are given in the appendix. 

STAGS model .- The single finite element type used in the STAGS model of this and the 
other examples is a six- node, triangular, combined membrane and bending element. The 
element is based on the Clough- Felippa triangle and has a displacement formulation. Midside 
nodes allow a single element across the depth of a blade stiffener to represent accurately its 
overall inplane bending behavior. The STAGS designation for this element is 422. 
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(b) Oblique view of model. 

Figure 12.- EAL finite element model for example 1 
composite blade-stiffened panel. 
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The finite element grid chosen for the STAGS model is shown in figure 13. The total 
number of triangular elements in the model is 2000. The total number of degrees of freedom 
in the model is 21 541. A more refined model containing four elements between blades was 
used as a check for selected loadings. 

Results .- Buckling results obtained with PASCO, EAL, and STAGS are shown in fig- 
ure 14. The curves indicate VIP ASA and smeared stiffener solutions, and the symbols indi- 
cate EAL and STAGS solutions. The solid curve represents the VIP ASA solution for buckling 
half-wavelength X equal to L. The dotted line at the top of the figure represents the 
VIP ASA solution for X equal to L/2. Rather than presenting VIP ASA solutions for only the 
lowest buckling load, VIP ASA solutions are presented for the X = L mode and for the lowest 
mode other than the X = L mode. In this example, the additional mode is the X = L/2 
mode. The X = L mode is always shown because its conservatism in the presence of shear 
is a focus of this report. ___ 

The dashed curve in figure 14 represents the smeared stiffener solution and indicates 
solutions for the lowest buckling load of the set X = W, W/2, W/3, ..., where W is the panel 
width. The slope change in the dashed curve that occurs at N x equal to approximately 
750 lb/in. (131 kN/m) indicates a change in mode shape for the smeared stiffener solution. 
For N x less than 750 lb/in., the buckling half-wavelength transverse to the stiffeners is 
equal to 15.0 in. (38.1 cm), which is 3 times the stiffener spacing (5.0 in.). For N x greater 
than 750 lb/in., the buckling half-wavelength transverse to the stiffeners is equal to 30 in., 
which is 6 times the stiffener spacing. 

In this example, the smeared stiffener solution gives reasonably accurate estimates of 
the solution for all combinations of N x and N^. For the pure shear loading, the smeared 
stiffener solution is about 5 percent lower than the EAL solution. For this same loading, the 
VEP ASA solution for X = L is about 63 percent lower than the EAL solution. For the pure 
compression loading, the VIP ASA solution for X = L and the EAL solution agree within 
0.3 percent. For the pure shear loading, the STAGS and EAL solutions agree within 
0.5 percent. 

Detailed comparisons and benchmark calculations for six loadings are presented in 
table IV. In this table, the quantity denoted FACTOR is the buckling load in terms of a scale 
factor for the specified loading. For example, for the loading N x = 2000 lb/in., 

Njjy = 1000 lb/in. (N x = 350.3 kN/m, Njjy = 175.1 kN/m), the EAL solution of 
FACTOR = 0.4764 means that the buckling load is N x = 0.4764 x 2000 = 952.8 lb/in. 

(166.9 kN/m), N xy = 0.4764 x 1000 = 476.4 lb/in. (83.43 kN/m). 

Also shown in table IV is an indication of the buckling mode shape determined by EAL 
for each loading. The designation of mode shapes as either symmetric (Sym) or antisym- 
metric (Anti) depends upon whether the buckling displacement w is symmetric or anti- 
symmetric with respect to a point at the center of the panel. 
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(a) Three views of model. 



(b) Oblique view of model. 

Figure 13.- STAGS finite element model for example 1 3 composite 

blade-stiffened panel. 
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Figure 14.- Buckling load interaction for example 1, 
composite blade-stiffened panel. 


Finally, the buckling mode shape obtained with EAL for the pure shear loading is shown 
in figure 15. The contour plot (fig. 15(b)) represents the buckling displacement w of the 
panel skin. This plot shows that the buckling half-wavelength transverse to the stiffeners is 
approximately equal to 3 times the stiffener spacing, as predicted by the smeared stiffener 
solution. This contour plot is an illustration of an antisymmetric buckling mode. The buck- 
ling mode shapes obtained with PASCO and EAL for pure compression are shown in 
figure 16. 
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(a) Mode shape obtained with EAL. 


Undeformed shape 

Buckling mode shape 



(b) Mode shape obtained with PASCO. 

Figure 16.- Buckling mode shape for pure longitudinal compression 
for example 1, composite blade-stiffened panel. 
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TABLE I.- WALL CONSTRUCTION FOR EACH PLATE ELEMENT IN 
EXAMPLE 1, COMPOSITE BLADE-STIFFENED PANEL 

Only half the laminate is defined for each plate elemenf 
because all laminates are symmetric _ 


Layer number, 
starting with 
outside layer 

Thickness 

Fiber 

orientation, 
6 , deg 

in. 

cm 


Plate elements 1 and 3 


1 

0.00550 

0.01397 

45 

2 

.00550 

.01397 

-45 

3 

.00550 

.01397 

-45 

4 

.00550 

.01397 

45 

5 

.00550 

.01397 

0 

6 

.04950 

.12573 

90 

Plate element 2 

1 

0.00550 

0.01397 

45 

2 

.00550 

.01397 

-45 

3 

.00550 

.01397 

-45 

4 

.00550 

.01397 

45 

5 

.01100 

.02794 

0 


TABLE H.- PROPERTIES OF MATERIALS USED IN CALCULATIONS 


Graphite- epoxy: 
Ej, psi (GPa) 
E 2 , psi (GPa) 
E 12 , psi (GPa) 

^1 

M2 

Aluminum: 

E, psi (GPa) 

M 


. 19 x 10 6 

(131) 

1.89 x 10 6 

(13.0) 

0.93 x 10 6 

(6.41) 


0.38 


0.0378 

10.5 x 10 6 

(72.4) 


0.32 


TABLE III.- PREBUCKLING LOAD DISTRIBUTION FOR EACH PLATE 
ELEMENT IN EXAMPLE 1, COMPOSITE 
BLADE-STIFFENED PANEL 

[Applied loading is N x = N^ = 1000.00 lb/in. (175.13 kN/m)] 


Load 

type 

Internal load distribution 

Plate elements 1 and 3 

Plate element 2 


lb/ in. 

kN/m 

lb/in. 

kN/m 

N x 

811.35 

142.09 

697.67 

122.18 

N xy 

1000.00 

175.13 

0 

0 


TABLE IV.- VALUES OF FACTOR TO OBTAIN BUCKLING LOADS FOR 
EXAMPLE 1, COMPOSITE BLADE-STIFFENED PANEL 


Applied load 

Factor 

EAL 

mode 

shape 

N x 

N xy 

VIPASA 

Smeared 

stiffener 

solution 

EAL 

STAGS 

lb/ in. 

kN/m 

lb/ in. 

kN/m 

11 

11 

< 

oo 

0 

0 

1000 

175.1 

0.5706 

1.6600 

1.4683 

1.5525 

1.5565 

Anti 

200 

35.0 

1000 

175.1 

.5339 

1.5577 

1.3098 

1.3985 


Anti 

500 

87.6 

1000 

175.1 

.4851 

1.4217 

1.1222 

1.2060 


Anti 

1000 

175.1 

1000 

175.1 

.4173 

1.2332 

.8222 

.8397 


Sym 

2000 

350.3 

1000 

175.1 

.3194 


.4690 

.4764 


Sym 

1000 

175.1 

0 

0 

1.0005 


.9970 

1.0030 


Sym 
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Example 2 - Metal Blade-Stiffened Panel 


Panel description .- A repeating element of the metal blade- stiffened panel is shown in 
figure 17. Element widths and thicknesses are also shown. The material properties used in 
the calculations are E = 10.5 x 10^ psi (72.4 GPa), p = 0.32. The prebuckling load distri- 
bution is given in table V. The orthotropic stiffnesses used for the smeared stiffener solution 
are = 89 074 lb- in., D 2 2 = 578.06 lb-in., and D33 = 314.61 lb- in. 



Figure 17.- Repeating element for example 2, metal blade-stiffened panel. 


***** EXAMPLE 2, METAL BLADE- STIFFENED PANEL ***** 

$CONDAT 

$ 

$PANEL 

B=2 .5 , 1.352, 2.5, 

T>=.042, .029, 

THET= 0, 0, 

KWALL( 1 , 1)= 1 , 

KWALL( 1 , 2)=2 , 

IWALL= 1,2,1, 

HCARJ>4,-4,2,90,0, 

2,121,4, 

4.5.1.3, -121, 

ICARI>=5 ,1,3,1,-909,0900, 

3. 2. 3. 4, 

3, 3, 4, 3, 

3.4, -909,0900, 

ICREB=6, 

N0BAY= 6 , 

SHEAR= 1 . , 

EL=30, 

MINLAM= 30 , 

NLAM= 1,2,3, 

IBC= 1 , 

IP=2, 

NX=1000. , 

NXY=-1000. , 

$ 

$MATER 

E1=10.5E6 , E2=10.5E6, E12=3.9772727E6 , ANU1=.32, RHO=.l, 

$ 

Figure 18.- Sample PASCO input for example 2, metal blade-stiffened panel. 
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PASCO input .- Sample PASCO input for this example is shown in figure 18. As is the 
case with all examples in this report, ICARD input is included to get detailed plots of buckling 
mode shapes. The same numerical results are obtained when the modeling is carried out 
with HCARD input only. The repeating element shown in figure 17 was generated with PASCO 
input. 



0 100 200 300 400 500 600 700 800 900 1000 


Figure 19.- Buckling load interaction for example 2, 
metal blade-stiffened panel. 
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EAL and STAGS models .- The same finite element type and grid pattern that are used in 
the EAL model of the first example are used in the EAL model for this example. Also, the 
finite element type and grid pattern used in the STAGS model of the first example are used in 
the STAGS model of this example. 

Results .- Buckling results are shown in figure 19. The same general approach used for 
presenting the buckling results in the first example is used in this and subsequent examples. 
Curves indicate PASCO results, and symbols indicate EAL and STAGS results. 

A portion of the standard VIP ASA solution for X = L/2, shown by the dotted line near 
the top of figure 19, is below the smeared stiffener solution. The slope changes in the 
smeared stiffener curve that occur at N x equal to approximately 730 lb/ in. (130 kN/ m) and 
950 lb/in. (166 kN/m) indicate a change in mode shape. For values of N x less than about 
730 lb/in., the buckling half-wavelength transverse to the stiffeners is equal to 10.0 in. 

(25.4 cm), which is twice the stiffener spacing. For values of N x greater than 730 lb/ in. 
but less than 950 lb/in., the buckling half- wavelength transverse to the stiffeners is 15.0 in. 
(38.1 cm), which is 3 times the stiffener spacing. For values of N x greater than 950 lb/in., 
the buckling half-wavelength transverse to the stiffeners is 30.0 in. (76.2 cm), which is 
6 times the stiffener spacing. 

The curves representing the VIP ASA results for X = L/2 and the smeared stiffener 
solution are above the EAL and STAGS results. For the pure shear case, an examination of 
the EAL buckling mode shape (fig. 20) shows that the buckling mode is an overall mode (X = L) 
rather than a X = L/2 mode, which might have been assumed because the X = L/2 solution 
is near the EAL and STAGS solutions. One possible factor contributing to the error in the 
smeared stiffener solution near N x = 0 is that the buckling half-wavelength transverse to 
the stiffeners is too short to be a valid solution. A buckling half-wavelength that is at least 
2.5 times the stiffener spacing would, in general, be necessary for a valid solution. However, 
even when such a buckle mode shape requirement is met, the smeared stiffener solution may 
not provide a good prediction of the buckling load. These results and results for subsequent 
examples show the danger of depending upon a smeared stiffener solution - even if the con- 
servative assumption is made that the panel width is infinite rather than finite. 

Detailed comparisons of solutions from PASCO, EAL, and STAGS for six loadings are 
presented in table VI. The buckling mode shape obtained with PASCO for the pure compres- 
sion loading is shown in figure 21. 
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Figure 20.- Buckling mode shape for pure shear for example 2> 
metal blade-stiffened panel. Obtained with EAL. 


— Undeformed shape 

Buckling mode shape 



Figure 21.- Buckling mode shape for pure longitudinal compression 
for example 2 5 metal blade-stiffened panel. Obtained with 
PASCO. 


TABLE V.- PREBUCKLING LOAD DISTRIBUTION FOR EACH PLATE 
ELEMENT IN EXAMPLE 2, METAL BLADE-STIFFENED PANEL 

^Applied loading is N x = N^ = 1000.00 lb/ in. (175.13 kN/m)J 




Internal load distribution 


Load 

type 

Plate elements 1 and 3 

Plate element 2 

lb/in. 

kN/m 

lb/in. 

kN/m 

N x 

842.67 

147.57 

581.84 

101.90 

N xy 

1000.00 

175.13 

0 

0 


TABLE VI.- VALUES OF FACTOR TO OBTAIN BUCKLING LOADS FOR 
EXAMPLE 2, METAL BLADE-STIFFENED PANEL 


Applied load 

Factor 

EAL 

mode 

shape 

N x 

£ 

VIP ASA 

Smeared 

stiffener 

solution 

EAL 

STAGS 

lb/ in. 

kN/m 

lb/ in. 

kN/ m 

« 

/< 

X = L/2 

0 

0 

1000 

175.1 

0.3118 

0.8450 

0.9823 

0.8138 

0.8179 

Sym 

400 

70.0 

1000 

175.1 

.2877 

.7742 

.8423 

.7195 


Sym 

1000 

175.1 

1000 

175.1 

.2568 

.6849 

.6879 

.6061 


Sym 

2000 

350.3 

1000 

175.1 

.2159 


.4637 

.4444 


Anti 

5000 

875.6 

1000 

175.1 

.1413 


.1975 

.1929 


Sym 

1000 

175.1 

0 

0 

.9710 

i 


.9969 

.9759 


Sym 
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Example 3 - Heavily Loaded, Composite Blade-Stiffened Panel 


Panel description .- A repeating element of the heavily loaded, composite blade- stiffened 
panel is shown in figure 22. The wall construction for each plate element is given in 
table VH. Only half the laminate is defined for each plate element because all laminates are 
symmetric. Plate element numbers are indicated by the circled numbers in figure 22. 



Figure 22.- Repeating element for example 3» heavily loaded, 
composite blade-stiffened panel. 


***** EXAMPLE 3, COMPOSITE BLADE-STIFFENED PANEL, HEAVILY LOADED ***** 
$CONDAT 
$ 

$PANEL 
GRANGER 10, 

B= 2 .5 , 1.97, 2.5, 

T=* .00637 , .0249, .0416, .00823, .0675, 

THET=45 , 0, 90, 45, 0, 

KWALL( 1,1)=1,-1,-1,1,2,3, 

KWALL( 1,2) =4, -4, -4, 4, 5, 

IWALW,2,1, 

HCARI>=4 ,-4,2,90,0, 

2,121,4, 

4.5.1.3, -121, 

ICARD=5, 1,3, 1,-909, 0900, 

3. 2. 3. 4, 

3,3 ,4,3, 

3.4, -909,0900, 

N0BAY=6 , 

ICREP= 6 , 

ELf30 , 

MINLAtt=30, 

NLAM= 1,2,3, 

SHEAR= 1 . , 

IBOl, 

11^=2, 

NX= 1000 . , 

NXY=- 1000 . , 

$ 

$MATER 

El= 19 .E6 , E2= 1 .89E6 , E12=.93E6, ANU1=.38, RHO=.0571, 

$ 
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Figure 23.- Sample PASCO input for example 3, heavily loaded 
composite blade-stiffened panel. 



Values for material properties are the same as those used in example 1 and are given in 
table n. The prebuckling load distribution is given in table VDI, and the orthotropic stiff- 
nesses are D-q = 958 850 lb-in., D 22 = 3339.4 lb-in., and D 33 = 4066.7 lb-in. 

PASCO input and EAL and STAGS models .- Sample PASCO input for this example is 
shown in figure 23. The same finite element types and grid patterns used in the EAL and 
STAGS models of the first example are used in this example. 

Results .- Buckling solutions for example 3 are shown in figure 24. The solid curve 
indicates the VIPASA solution for X = L. The dotted curves indicate VIPASA solutions for 



Figure 24.- Buckling load interaction for example 3, heavily loaded 

composite blade-stiffened panel. 
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X = L/2, L/4, and L/5. The dashed curve represents the smeared stiffener solution. As in 
the first example, the slope changes in the dashed curve indicate changes in mode shape. 
These changes in slope occur at N x equal to approximately 4000 lb/ in. (700 kN/m), 

9200 lb/in. (1610 kN/m), and 10 300 lb/in. (1800 kN/m). For N x less than about 
4000 lb/in., the buckling half-wavelength transverse to the stiffeners is 1.5 times the 
stiffener spacing. For N x greater than about 4000 lb/in. but less than about 9200 lb/in., 
the buckling half-wavelength transverse to the stiffeners is 2 times the stiffener spacing. 

For N x greater than about 9200 lb/ in. but less than about 10 300 lb/ in., the buckling half- 
wavelength transverse to the stiffeners is 3 times the stiffener spacing. 

The curves representing the VIP ASA results for X = L/2, L/4, and L/5 and the 
smeared stiffener solution are all above the EAL and STAGS results. As in example 2, the 
buckling mode for the EAL analysis with a pure shear loading is an overall mode - not a 
X = L/2 mode. 

Detailed comparisons of solutions from PASCO, EAL, and STAGS for six loadings are 
presented in table IX. The buckling mode shape obtained with EAL for the pure shear loading 
is shown in figure 25. The buckling mode shape obtained with PASCO for the pure compres- 
sion loading is shown in figure 26. 



Figure 25.- Buckling mode shape for pure shear for example 3 a 
heavily loaded, composite blade-stiffened panel. Obtained 
with EAL. 
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Undeformed shape 

Buckling mode shape 



Figure 26.- Buckling mode shape for pure longitudinal compression 
for example 3, heavily loaded, composite blade-stiffened panel. 
Obtained with PASCO. 


TABLE VH.- WALL CONSTRUCTION FOR EACH PLATE ELEMENT IN 
EXAMPLE 3, HEAVILY LOADED, COMPOSITE 
BLADE-STIFFENED PANEL 

Only half the laminate is defined for each plate element 
because all laminates are symmetric 


Layer number, 
starting with 
outside layer 

Thickness 

Fiber 

orientation, 
6 , deg 

in. 

cm 


Plate elements 1 and 3 


1 

0.00637 

0.01618 

45 

2 

.00637 

.01618 

-45 

3 

.00637 

.01618 

-45 

4 

.00637 

.01618 

45 

5 

.02490 

.06325 

0 

6 

.04160 

.10566 

90 

Plate element 2 

1 

0.00823 

0.02090 

45 

2 

.00823 

.02090 

-45 

3 

.00823 

.02090 

-45 

4 

.00823 

.02090 

45 

5 

.06750 

.17145 

0 



TABLE m- PREBUCKLING LOAD DISTRIBUTION FOR EACH 
PLATE ELEMENT IN EXAMPLE 3, HEAVILY LOADED, 
COMPOSITE BLADE-STIFFENED PANEL 

^Applied loading is N x = N^ = 1000.00 lb/in. (175.13 kN/m)J 



Internal load distribution 

Load 

type 

Plate elements 1 and 3 

Plate element 2 

lb/in. 

kN/m 

lb/ in. 

kN/m 

N x 

559.71 

98.02 

1117.50 

195.70 

N xy 

1000.00 

175.13 

0 

0 


TABLE IX.- VALUES OF BUCKLING LOADS FOR EXAMPLE 3, HEAVILY 
LOADED, COMPOSITE BLADE-STIFFENED PANEL 


Applied load 



Factor 



EAL 

mode 

shape 

N 

X 

N xy 

VIPASA 

Smeared 

stiffener 

solution 

EAL 

STAGS 

lb/ in. 

kN/m 

lb/ in. 

kN/m 

X = L 

X - L/2 

0 

0 

1000 

175.1 

2.9179 

6.6855 

9.2435 

6.4424 

6.470 

Anti 

500 

87.6 

1000 

175.1 

2.6703 

6.0262 

8.0628 

5.753 


Anti 

1000 

175.1 

1000 

175.1 

2.4540 

5.4549 

6.7945 

5.1630 


Anti 

2000 

350.3 

1000 

175.1 

2.0972 

4.5289 

4.8627 

4.124 


Sym 

4000 

700.5 

1000 

175.1 

1.5949 


2.6424 

2.4543 


Sym 

1000 

175.1 

0 

0 

9.9724 


10.7300 

10.076 


Sym 


32 



Example 4 - Metal Blade -Stiffened Panel With Thin Skin 

Panel description .- A repeating element of the metal blade- stiffened panel with thin skin 
is shown in figure 27. Element widths and thicknesses are also shown. The material prop- 
erties used in the calculations are E = 10.5 x 10^ psi (72.4 GPa), p = 0.32. Except for the 
thickness of the skin, this panel is the same as that used in example 2. The prebuckling load 
distribution is given in table X. The orthotropic stiffnesses are = 82 490 lb-in., 

D22 = 121.92 lb-in., and D33 = 86.641 lb-in. 



Figure 27*- Repeating element for example 4, metal 
blade-stiffened panel with thin skin. 

PASCO input and EAL and STAGS models .- PASCO input for this example is the same as 
that given in figure 18 except that thickness T(l) = 0.025. The same finite element types and 
grid patterns used in the EAL and STAGS models of the first example are used in this 
example. 

Results .- Buckling results obtained with PASCO, EAL, and STAGS are shown in fig- 
ure 28. The same general approach used for presenting buckling data in previous examples 
is used in this example. Because the skin is thinner for this example than for example 2, the 
dominant buckling mode is a local mode rather than an overall mode. 

The solid curve represents the standard VIP ASA solution for X = L. The dotted curve 
represents the standard VIPASA solutions for X = L/5 and X = L/6. The dashed curve 
represents the smeared stiffener solution. For the pure shear loading, the smeared stiffener 
solution has a buckling half-wavelength transverse to the stiffeners that is 1.2 times the 
stiffener spacing. That buckling half-wavelength is much too short to provide a valid 
solution. 
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Figure 28.- Buckling load interaction for example 4, metal 
blade-stiffened panel with thin skin. 


The EAL results are very near the VIP ASA short-wavelength results indicated by the 
dotted curve. Except for the pure shear case, the EAL results are slightly higher than the 
VIPASA short-wavelength results. For the pure shear case, the lowest eigenvalue for both 
EAL and STAGS appears to be primarily an overall mode (in the stiffener direction) rather 
than a local mode. The mode shape for this first eigenvalue is shown in figure 29. The 
second eigenvalue obtained from EAL for a pure shear loading is a local mode and is shown 
in figure 30. 

Detailed comparisons of solutions from PASCO, EAL, and STAGS for six loadings are 
presented in table XI. The buckling mode shapes obtained with EAL and PASCO for the pure 
compression loading are shown in figure 31. 
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Figure 30.- Buckling mode shape for second eigenvalue, 
pure shear, for example 4, metal blade-stiffened 
panel with thin skin. Obtained with EAL. 


Undeformed shape 

Buckling mode shape 



(a) Mode shape obtained with PASCO. 

Figure 31.- Buckling mode shape for pure longitudinal 
compression for example 4, metal blade-stiffened 
panel with thin skin. 




Because the first eigenvalue for pure compression and the second eigenvalue for pure 
shear are local modes, a simplified local buckling analysis was also performed on this panel. 
For the pure compression case, local buckling of the skin can be estimated by assuming that 
the skin between stiffeners is a narrow, simply supported strip. The buckling load of the 
skin is given by 



(3) 


where D is the plate bending stiffness (121.9 lb-in.) and b is the plate width (5 in.). The 
result is = 192 lb/in. According to PASCO, when the panel buckles, the load in the skin 

is 225 lb/in. In this case, the stiffeners are providing rotational support to the skin. This 
can be shown by calculating the buckling load of the stiffeners with a formula (ref. 16, 
pp. 360-362) that is similar to equation (3). That formula is 




(4) 


where D is 190.2 lb-in. and b is 1.352 in. The result is N v = 513 lb/in. According 

A cr 

to PASCO, when the panel buckles, the load in the stiffeners is 261 lb/in. The skin is destabi- 
lizing the stiffeners, or, conversely, the stiffeners are stabilizing the skin. Further discus- 
sion of the concept of rotational restraint at the juncture of plate elements is presented in 
reference 17. (See, particularly, fig. 14 of ref. 17.) 

For the pure shear case, local buckling of the skin can be approximated by the formula 
(ref. 16, p. 383) 


N- 


xy 


= 5.35 


cr 


tt 2 D 

b 2 


(5) 


where D and b are given following equation (3). The result is N—., = 257 lb/in. 

xy cr 

According to PASCO and EAL, the shear load in the skin when the panel buckles in a local 
mode is about 300 lb/ in. As in the pure compression case, the stiffeners are providing 
rotational support to the skin. 
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TABLE X.- PREBUCKLING LOAD DISTRIBUTION FOR EACH PLATE 
ELEMENT IN EXAMPLE 4, METAL BLADE -STIFFENED 
PANEL WITH THIN SKIN 

[Applied loading is N x = N X y = 1000.00 lb/ in. (175.13 kN/m)J 




Internal load distribution 


Load 

type 

Plate elements 1 and 3 

Plate element 2 

lb/in. 

kN/m 

lb/in. 

kN/m 

N x 

761.23 

133.31 

883.03 

154.64 

N xy 

1000.00 

175.13 

0 

0 


TABLE XI.- VALUES OF FACTOR TO OBTAIN BUCKLING LOADS 
FOR EXAMPLE 4, METAL BLADE-STIFFENED 
PANEL WITH THIN SKIN 


Applied load 

Factor 

EAL 

mode 

shape 

N 

X 

N xy 

VIP ASA 

Smeared 

stiffener 

solution 

EAL 

STAGS 

lb/ in. 

kN/ m 

lb/ in. 

kN/ m 

X = L 

X = L/5, L/6 

0 

0 

1000 

175.1 

0.1761 

0.2961 

0.6062 

0.2767 

0.2773 

Anti 

0 

0 

1000 

175.1 




.305 a 


Sym 

400 

70.0 

1000 

175.1 

.1671 

.2428 

.5541 

.2491 


Anti 

1000 

175.1 

1000 

175.1 

.1548 

.1840 

.4760 

.1881 


Anti 

2000 

350.3 

1000 

175.1 

.1374 

.1227 

.3695 

.1253 


Anti 

4000 

700.5 

1000 

175.1 

.1112 

.06984 

.2222 

.07064 


Anti 

1000 

175.1 

0 

0 

.8611 

.2958 

.9097 

.2965 


Anti 


a Second eigenvalue. 


39 



Example 5 - Composite Hat-Stiffened Panel 

Panel description .- A repeating element of the composite hat-stiffened panel is shown 
in figure 32. Plate element widths are also shown. The wall construction for each plate ele- 
ment is given in table XH. Only half the laminate is defined for each plate element because 
all laminates are symmetric. Plate element numbers are indicated by the circled numbers 
in figure 32. Values for material properties are the same as those used in example 1 and are 
given in table II. The prebuckling load distribution is given in table XHI. The orthotropic 
stiffnesses are D-j^ = 242 270 lb-in., D22 = 447.77 lb-in., and D33 = 8511.2 lb-in. 



Figure 32.- Repeating element for example 5, 
composite hat-stiffened panel. 

PASCO input and EAL model .- Sample PASCO input for this example is shown in fig- 
ure 33. (Note that the PASCO input parameter SHEAR is set to 0.3, which indicates that for 
the smeared stiffener solution, the twisting stiffness D33 has been multiplied by 0.3. The 
reduction is used because the formula for D33 is based on the enclosed area of the closed 
cross section, which, by experience, tends to overestimate the effective stiffness. The value 
for D33 given above is the value after multiplying by 0.3. See refs. 11 and 12.) The finite 
element type used in the EAL model is the same as that used in the previous models. The 
finite element grid chosen for the EAL model is shown in figure 34. There are 36 elements 
along the length. One finite element is used for plate elements 1, 3, and 5, and two finite 
elements are used for plate elements 2, 4, and 6. There are 1944 finite elements and 
1813 nodes. 
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***** 


EXAMPLE 5, COMPOSITE HAT- STIFFENED PANEL 
$ CONDAT 
$ 

$PANEL 
GRANGER 10, 

EIf30 , 

B=1 .3 , 2.4, 1.3, 1.3, 1.1, 1.3, 

T= .010315, .009953, .016955, .025383, 

THET=45 , 0, 0, 0, 

KWALL( 1 , 1)= 1,-1, -1,1, 2, 

KWALL( 1,2)= 1,-1, 3, 

KWALL(1,3)= 1,-1, 

KWALL( 1,4)= 1,-1, 4, 

IWALL= 1,2, 1,3, 4, 3, 

HCARD=4,-16 ,6,60,0, 

4 ,-14,4 ,-60,0, 

4.7.14.5.16, 

5. 8.1. 2, -7, 3, 

ICARI>=5, 1,2, 1,-909, 0900, 

3. 2. 5. 2, 

3,2,3,14, 

3, 3, 4, 5, 

3.4.5.16, 

3. 5. 6. 3, 

3,6,-909,0900, 

ICREP=6 , NOBAY= 6 , 

MINLAM=30, 

NLAM= 1,2, 

SHEAR= . 3 , 

IBC=1, 

NX=1000. , 

NXY=-1000 . , 

$ 

$MATER 

E1=19.E6, E2=1.89E6, E12=.93E6, ANU1=.38, 

$ 


***** 


RHO.0571, 


Figure 33.- Sample PASCO input for example 5 
composite hat-stiffened panel. 






VIPASA, X = L 
VIPASA, X = LIZ 
Smeared stiffener solution 





Stiffener direction 






The curve representing the smeared stiffener solution is above the EAL results for the 
entire range of loadings. If D 33 had not been reduced by the factor 0.3, the discrepancy 
would have been even greater. The EAL result for pure shear is an overall mode. 

Detailed comparisons of solutions from PASCO and EAL for six loadings are presented 
in table XIV. The buckling mode shape obtained with EAL for the pure shear loading is shown 
in figure 36. The buckling mode shapes obtained with EAL and PASCO for the pure compres- 
sion loading are shown in figure 37. 


Undeformed shape 

Buckling mode shape 



(a) Mode shape obtained with PASCO. 



(b) Oblique view of mode shape obtained with EAL. 

Figure 37 . - Buckling mode shape for pure longitudinal- 
compression for example 5, composite hat-stiffened 
panel . 
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TABLE XH.- WALL CONSTRUCTION FOR EACH PLATE ELEMENT IN 
EXAMPLE 5, COMPOSITE HAT-STIFFENED PANEL 

Only half the laminate is defined for each plate element 
because all laminates are symmetric 


Layer number, 
starting with 
outside layer 

Thickness 

Fiber 

orientation, 
9, deg 

in. 

cm 


Plate elements 1 and 3 


1 

0.010315 

0.026200 

45 

2 

.010315 

.026200 

-45 

3 

.010315 

.026200 

-45 

4 

.010315 

.026200 

45 

5 

.009953 

.025281 

0 


Plate element 2 


1 

0.010315 

0.026200 

45 

2 

.010315 

.026200 

-45 

3 

.016955 

.043066 

0 

Plate elements 4 and 6 

1 

0.010315 

0.026200 

45 

2 

.010315 

.026200 



-45 

Plate element 5 

1 

0.010315 

0.026200 

45 

2 

.010315 

.026200 

-45 

3 

.025383 

.064473 

0 


TABLE Xm.- PREBUCKLING LOAD DISTRIBUTION FOR EACH PLATE 
ELEMENT IN EXAMPLE 5, COMPOSITE HAT-STIFFENED PANEL 

[Applied loading is N x = N X y = 1000.00 lb/ in. (175.13 kN/m)] 



Internal load distribution 

Load 

type 

Plate elements 1 and 3 

Plate element 2 

Plate elements 4 and 6 

Plate element 5 

lb/in. 

kN/m 

lb/ in. 

kN/m 

lb/ in. 

kN/m 

lb/ in. 

kN/m 

N x 

631.51 

110.59 

764.63 

133.91 

128.89 

22.57 

1079.87 

189.11 

N xy 

1000.00 

175.13 

626.98 

109.80 

373.02 

65.33 

373.02 

L . . 

65.33 


TABLE XIV.- VALUES OF FACTOR TO OBTAIN BUCKLING LOADS 
FOR EXAMPLE 5, COMPOSITE HAT-STIFFENED PANEL 


Applied load 

Factor 

EAL 

mode 

shape 

N x 

N xy 

VIPASA 

Smeared 

stiffener 

solution 

EAL 

lb/in. 

kN/m 

lb/in. 

kN/m 

A = L 

A = L/2 

0 

0 

1000 

175.1 

1.2906 

3.5433 

3.5698 

3.192 

Sym 

300 

52.5 

1000 

175.1 

1.2045 

3.3203 

3.1933 

2.932 

Sym 

600 

105.1 

1000 

175.1 

1.1258 

3.1147 

2.8822 

2.680 

Anti 

1000 

175.1 

1000 

175.1 

1.0313 

2.8656 

2.4229 

2.3268 

Sym 

2000 

350.3 

1000 

175.1 

.8385 

2.3525 

1.4265 

1.4062 

Sym 

1000 

175.1 

0 

0 

2.9952 


3.0351 

3.0042 

Sym 
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Example 6 - Composite Corrugated Panel 

Panel description .- A repeating element of the composite corrugated panel is shown in 
figure 38. Element widths are also shown. The wall construction for each plate element is 
given in table XV. Plate element numbers are indicated by the circled numbers in figure 38. 
Only half the laminate is defined for each plate element because all laminates are symmetric. 
Values for material properties are the same as those used in example 1 and are given in 
table n. The prebuckling load distribution is given in table XVI. The orthotropic stiffnesses 
are = 130 930 lb-in., D22 = 216.53 lb-in., and D33 = 184.74 lb-in. 



Figure 38 .- Repeating element for example 6 , composite corrugated panel. 


***** EXAMPLE 6, COMPOSITE CORRUGATED PANEL ***** 

$ CON DAT 

$ 

$PANEL 
ELr=30 , 

B= .73223, 1.46446, 1.46446, 1.46446, .73223, 

T=. 005479, .016836, 

THET-45 , 0, 

KWALL( 1,1)= 1,-1, -1,1, 2, 

KWALL(1,2)= 1,-1,-1,1, 

IWALL- 1,2, 1,2,1, 

HCARI>4 ,-6 , 2,-450 , -1, 

4,-7,4,450, -1, 

6. 8. 1.6. 3. 7. 5, 

ICARI>x 5 , 1 , 2 , 1 , -909 , 0900 , 

3. 2. 3. 6, 

3, 3, 4, 3, 

3. 4. 5. 7, 

3. 5. 6. 5, 

3,6,-909,0900, 

ICREP-6 , NOBAY- 6, 

MINLAM-30 , 

NLAM- 1,2, 

SHEAR- 1 . , 

IBOl, 

IP-2, 

NX- 1000., NXY— 1000. , 

$ 

$MATER 

E1-19.E6, E2-1.89E6, E12-.93E6, ANU1-.38, RH0-.0571, 


Figure 39*- Sample PASCO input for example 6, composite corrugated panel. 
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P ASCO input an d E AL model . - Sample PASCO input for this example is shown in fig- 
ure 39. The same finite element type used in the EAL model of the first example is used in 
this example. For the EAL finite element model, shown in figure 40, 36 elements are used 
down the length of the panel, and 2 elements are used across the width of each plate element. 




(a) Three views of model. 
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(b) Oblique view of model. 

Figure llO.— EAL finite element model for example 6, 
composite corrugated panel. 


Results.- Buckling results obtained with PASCO and EAL for this example are shown in 
figure 41. The curve for the smeared stiffener solution falls below the EAL results; there- 
fore, for this case, the smeared stiffener approach is conservative. This conservatism exists 



Figure 4l.- Buckling load interaction for example 6 , 
composite corrugated panel. 


even though the smeared stiffener solution predicts that the buckling half-wavelength trans- 
verse to the stiffeners is relatively small. For example, for the pure shear case, the smeared 
stiffener solution predicts that the buckling half-wavelength transverse to the stiffeners is 
only 1.2 times the stiffener spacing, which, for this case, is the corrugation pitch. Results 
from the other examples suggest that a buckling half-wavelength of only 1.2 times the stiffener 
spacing would cause the smeared stiffener solution to be very unconservative; that is, the 
solution would overestimate the buckling load. The buckling mode shape for the pure shear 
case, figure 42, shows that buckling is limited to the area near the side edges. 
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Figure 42.- Buckling mode shape for pure shear for example 6, 
composite corrugated panel. Obtained with EAL. 

Because these results do not fit the pattern of results from the other examples, several 
additional calculations were made. First, the correct simple support boundary conditions 
were applied to the sides of the smeared stiffener model. However, imposing these boundary 
conditions increased the buckling load by only about 3 percent for the pure shear case. 

Next, it was thought that the location of the side supports on the EAL model was causing 
the EAL buckling load to be larger than it would be for supports at other locations. To 
explore that possibility, a second EAL model was constructed. That model is identical to the 
original model except that the supports are shifted one quarter period to the right, as shown in 
figure 43. In the new model, the supports are at the midheight of the cross section. Calcula- 
tions with the new model showed that, instead of being smaller, the buckling load for the pure 





Figure 43.- Corrugated panel modeled with supports 
at midheight of cross section. 


shear loading increased by about 15 percent. As can be seen in figure 44, the buckling mode 
shape for the second model extends over the entire panel rather than being limited to the area 
near the side edges, which was the case for the original model (fig. 42). 

Another possible reason the smeared stiffener solution underestimates the buckling load 
for the corrugated panel is that the orthotropic stiffnesses are underestimated. In particular, 
neglecting restraint of warping may be a factor in underestimating the twisting stiffness. 
However, the twisting stiffness formula used in the PASCO smeared stiffener solution 
(eq. (13), ref. 11) includes extra terms which make it larger than the twisting stiffness calcu- 
lated by many standard approaches, such as the formula on page 367 of reference 18. 
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Figure 44.- Buckling mode shape for pure shear for corrugated panel 
supported as shown in figure 43. Obtained with EAL. 

It can only be concluded that although a smeared stiffener solution provides an estimate 
of the buckling load for an overall mode, it sometimes overestimates the buckling load and 
sometimes underestimates the load. This simplified approach cannot be counted on to be 
consistent. 

Detailed comparisons of results from PASCO and EAL for six loadings are presented in 
table XVII. The buckling mode shapes obtained with EAL and PASCO for the pure compres- 
sion case are shown in figure 45. 



(a) Mode shape obtained with EAL. 

Undeformed shape 

Buckling mode shape 



(b) Mode shape obtained with PASCO. 

Figure 45.- Buckling mode shape for pure longitudinal compression 
for example 6, composite corrugated panel. 
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TABLE XV.- WALL CONSTRUCTION FOR EACH PLATE ELEMENT 
IN EXAMPLE 6, COMPOSITE CORRUGATED PANEL 

[ Only half the laminate is defined for each plate element 
because all laminates are symmetric 


Layer number, 
starting with 
outside layer 

Thickness 

Fiber 

in. 

cm 

orientation, 
6, deg 


Plate elements 1, 3, and 5 


1 

0.005479 

0.013917 

45 

2 

.005479 

.013917 

-45 

3 

.005479 

.013917 

-45 

4 

.005479 

.013917 

45 

5 

.016836 

.042763 

0 


Plate elements 2 and 4 


1 

0.005479 

0.013917 

45 

2 

.005479 

.013917 

-45 

3 

.005479 

.013917 

-45 

4 

.005479 

.013917 

45 



TABLE XVI.- PREBUCKLING LOAD DISTRIBUTION FOR EACH PLATE 
ELEMENT IN EXAMPLE 6, COMPOSITE CORRUGATED PANEL 

[Applied loading is N x = N xy = 1000.00 lb/ in. (175.13 kN/m)] 



Internal load distribution 

Load 

type 

Plate elements 1, 3, and 5 

Plate elements 2 and 4 

lb/in. 

kN/m 

lb/in. 

kN/m 

N x 

1448.89 

253.74 

258.22 

45.22 

N xy 

1000.00 

175.13 

1000.00 

175.13 


TABLE XVH.- VALUES OF FACTOR TO OBTAIN BUCKLING LOADS 
FOR EXAMPLE 6, COMPOSITE CORRUGATED PANEL 


Applied load 

Factor 

N x 

N xy 

VIPASA 

Smeared 

stiffener 

solution 

EAL 

lb/in. 

kN/m 

lb/ in. 

kN/m 

X = L 

0 

0 

1000 

175.1 

0.4371 

1.0011 

1.248 

500 

87.6 

1000 

175.1 

.4010 

.8894 

1.1395 

1000 

175.1 

1000 

175.1 

.3692 

.7780 

1.0223 

2000 

350.3 

1000 

175.1 

.3161 

.5971 

.7077 

5000 

875.6 

1000 

175.1 

.2131 

.2865 

.2958 

1000 

175.1 

0 

0 

1.4897 

1.4462 

1.4918 




Example 7 - Metal Z- Stiffened Panel 


Panel description .- A repeating element of the metal Z- stiffened panel is shown in fig- 
ure 46. Element widths and thicknesses are also shown. The material properties used in the 



Figure 46.- Repeating element for example 7* metal Z-stiffened panel. 


***** EXAMPLE 7, METAL Z-STIFFENED PANEL ***** 

$CONDAT 

$ 

$ PANEL 

B=2 .5 , 1.352, 2.5, .5, 

T= *042, .029, 

THET= 0, 0, 

KWALL( 1 , 1)= 1 , 

KWALL( 1 , 2) = 2 , 

IWALL= 1 , 2 , 1 , 2 , 

HCARD=4 ,-5,2,90,0, 

3.121.4.5, 

4.6.1.3, -121, 

ICARD=5 , 1,4,1 ,-909 ,0900, 

3 > 3 , 4 , 5 , 

3. 2.3. 4, 

3,4,5, 3, 

3.5, -909,0900, 

ICREP=6, 

N0BAY=6 , 

Eb=30, 

MINLAM=30 , 

NLAM= 1 , 2 , 5 , 6 , 

SHEAR= 1 • , 

IBC= 1 , 

IP=2, 

NX= 1000 . , 

NXY=-1000 . , 

$ 

$ MATER 

El= 10.5E6 , E2= 10 . 5E6 , E12=3 .9772727E6 , ANU1=.32, RHO=.l, 

$ 

$ 

Figure 47-*- Sample PASCO input for example 7 5 metal Z-stiffened panel. 
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calculations are E = 10.5 x 10® psi (72.4 GPa), p = 0.32. Except for the addition of the 
flange on the stiffener, this panel is the same as that used in example 2. The prebuckling 
load distribution is given in table XVTU. The orthotropic stiffnesses are 
Du = l" 7 ® 3S)C) lb-in., D22 = 578.07 lb-in., and D 33 = 324.11 lb-in. 

PASCO input and EAL model .- PASCO input for this example is presented in figure 47. 
The same basic finite element grid pattern used in the EAL model of the first example is 
used in this example. The flange (plate element 4) is modeled with a single element across 
the width. 

Results .- Buckling results obtained with PASCO and EAL are shown in figure 48. The 
same general approach used for presenting buckling data in previous examples is used in this 
example. Curves indicate PASCO results, and symbols indicate EAL results. The results 
for this example are similar, in many ways, to the results for example 4. 



Figure 48 .- Buckling load interaction for example 7 , metal Z-stiffened panel. 
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The solid curve represents the standard VIP ASA solution for A = L. The dotted curve 
represents the standard VIPASA solutions for X = L/2, L/6, and L/7. The dashed curve 
represents the smeared stiffener solution. For the pure shear loading, the smeared stiffener 
solution has a buckling half-wavelength transverse to the stiffeners that is 1.5 times the 
stiffener spacing. 

Detailed comparisons of solutions from PASCO and EAL for six loadings are presented 
in table XIX. The buckling mode shape obtained with EAL for the pure shear loading is shown 
in figure 49. Note that because of the flange, the stiffener is unsymmetric; therefore, 



(a) Oblique view. 
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(b) Contour plot. 

Figure 49-- Buckling mode shape for pure shear for example 7> 
metal Z-stiffened panel. Obtained with EAL. 



buckling mode shapes are neither perfectly symmetric nor perfectly antisymmetric. The 
buckling mode shapes obtained with PASCO and EAL for pure compression are shown in 
figure 50. 



(a) Mode shape obtained with PASCO. 



(b) Oblique view of mode shape obtained with EAL. 

Figure 50.- Buckling mode shape for pure longitudinal 
compression for example 7, metal Z-stiffened 
panel. 


59 






TABLE XVHI.- PREBUCKLING LOAD DISTRIBUTION FOR EACH PLATE 
ELEMENT IN EXAMPLE 7, METAL Z -STIFFENED PANEL 

[Applied loading is N x = N^ = 1000.00 lb/ in. (175.13 kN/m)] 



Internal load distribution 

Load 

type 

Plate elements 1 and 2 

Plate elements 3 and 4 

lb/in. 

kN/m 

lb/ in. 

kN/ m 

N x 

796.34 

139.46 

549.85 

96.29 

N xy 

1000.00 

175.13 

0 

0 


TABLE XIX.- VALUES OF FACTOR TO OBTAIN BUCKLING LOADS 
FOR EXAMPLE 7, METAL Z -STIFFENED PANEL 


Applied load 

Factor 

EAL 

mode 

shape 

N x 

N xy 

VIPASA 

Smeared 

stiffener 

solution 

EAL 

lb/ in. 

kN/m 

lb/in. 

kN/ m 

X = L 

X = L/2 

X = L/6, L/7 

0 

0 

1000 

175.1 

0.4924 

1.1282 


1.6050 

1.042 

Anti 

400 

70.0 

1000 

175.1 

.4606 

1.0336 


1.4343 

.948 

Anti 

1000 

175.1 

1000 

175.1 

.4186 

.9123 

0.8420 

1.1966 

.825 

Anti 

1000 

175.1 

500 

87.6 

.7224 


1.1292 

1.7501 

1.146 

Anti 

1000 

175.1 

200 

35.0 

1.2473 


1.3039 

1.9579 

1.315 

Sym 

1000 

175.1 

0 

0 

1.8662 


1.3503 

1.9768 

1.356 

Sym 
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RECENT DEVELOPMENTS 


During the final stages of preparation of this paper, a closely related paper dealing with 
the buckling analysis of shear-loaded stiffened panels was published. In that paper (ref. 19), 
a method is presented for enforcing desired support conditions along the ends of stiffened 
panels loaded in shear. The method combines Lagrangian multipliers, which are used to 
impose the support conditions, and the computer program VTPASA (refs. 1 to 3). The method 
described in reference 19 has been incorporated in a computer program denoted VICON. 

VICON retains many of the advantages and computational efficiencies of the VIP ASA 
program. VICON can calculate the buckling load of a stiffened panel with rather general rigid 
and elastic supports. In particular, VICON can be used when the buckling mode shape is 
required to have zero deflection along a rectangular boundary. In one calculation presented 
in reference 19, VICON is shown to give results in good agreement with EAL for example 2 of 
this paper. 


CONCLUDING REMARKS 

Buckling results obtained with several buckling analysis procedures are presented in 
this paper for seven stiffened panels. Four of the panels are blade- stiffened panels. One 
panel has a small flange at the tip of each blade to produce a Z- stiffened panel. One panel has 
hat stiffeners, and one panel is a corrugated panel. Four of the panels are made of graphite- 
epoxy; three panels are made of aluminum. The analysis procedures used to calculate buck- 
ling loads are the computer program PASCO, which includes both a linked-plate analysis pro- 
cedure denoted VIP ASA and a smeared stiffener solution, and two finite element computer 
programs, EAL and STAGS. For one panel, additional solutions for local buckling are 
obtained from simple formulas. 

The paper has two primary objectives: to examine the validity of the smeared stiffener 
solution and to provide a set of accurate benchmark calculations for evaluating other analysis 
procedures. A third objective is to provide insight into the buckling characteristics of 
stiffened panels. Numerous figures showing buckling mode shapes are included to help 
accomplish this third objective. 

A smeared stiffener solution is obtained by averaging or "smearing" the stiffnesses of 
the stiffeners over the panel. The result is an equivalent orthotropic plate that can be 
analyzed with relatively simple procedures. A smeared stiffener solution is only used when 
it is assumed that the buckling mode shape has a buckling wave that extends over several 
stiffeners. The buckling results presented in this paper indicate that the smeared stiffener 
solution should be used only with caution. Inaccuracies arise for at least two reasons. First, 
it may be difficult to calculate the proper orthotropic stiffnesses - particularly the twisting 
stiffness. Second, local stiffener deformations are ignored. These local deformations can 


62 



be important even when the buckling wave extends over several stiffeners. In the results 
presented in this paper, the smeared stiffener solution generally predicted a buckling load 
that was higher than the buckling load predicted by the more accurate analyses. 


Another conclusion of this study is that a PASCO solution can be very helpful when 
developing an adequate finite element model. Because PASCO is essentially exact for certain 
loadings and boundary conditions, the adequacy of a finite element solution (e.g., modeling 
detail, element capability, and eigenvalue solution) can be checked by comparing the finite 
element solution with the PASCO solution for a loading system and set of boundary conditions 
known to be exact in PASCO. Once the finite element procedure is determined to be adequate 
for the check case, the finite element analysis can be carried out with the desired loading and 
boundary conditions. Although this approach does not guarantee an accurate solution for the 
new loadings and boundary conditions, it can quickly point out deficiencies. 

Finally, these results are for buckling of perfect panels. All real panels have imper- 
fections that reduce their buckling loads. Analysis procedures can generally account for 
some of these imperfections. Also, buckling does not necessarily mean failure. For some 
panels and loadings, material strength failures occur before the buckling load is reached. 
Other panels may have substantial postbuckling strength. Still, an accurate prediction of the 
buckling load can be important in understanding the behavior of stiffened panels. 

Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
October 24, 1983 
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APPENDIX 


EAL RUNSTREAMS USED FOR PANEL ANALYSES 

The EAL runstreams used to carry out the panel studies are included in this section. 

An EAL runstream is somewhat analogous to a FORTRAN subroutine in that it consists of a 
collection of stored EAL input statements that can be processed by invoking a single command. 
Several of the runstreams are quite general and can be used in the analysis of any prismatic, 
plate-type structure composed of repeating elements. These runstreams were used for all 
panel configurations considered in this paper. To illustrate the use of these more general 
runstreams in performing a panel buckling analysis, the specific input is presented for 
example 1, the composite blade- stiffened panel. 

The input for the panel buckling analyses consists of a set of commands which initiate 
the execution of various EAL runstreams which, in turn, call other lower level runstreams. 
The execution flow and the runstreams for example 1 are described by starting with the high- 
est level and proceeding to the lower level, general purpose runstreams. 

The following commands control the execution of the example 1 buckling analysis: 

♦CALL (PANEL, BLADE) 

♦CALL (PANEL, GEOM) 

♦CALL (SA, BLADE) 

♦CALL (PANEL, BUCK) 

where the *CALL command merely begins execution of the runstream named between 
parentheses. 


Runstream PANEL, BLADE 

The runstream PANEL, BLADE, which is shown in figure Al, applies specifically to the 
particular problem being solved - in this case, the composite blade- stiffened panel (example 1) 
under a pure compression load. The main purpose of the runstream is to define several EAL 
tables used by the general purpose runstreams, which generate the panel geometry and directly 
apply the prebuckling stress state. The register NMNX defines the number of nodes to be 
generated along the panel length. The register NREP defines the number of repeating ele- 
ments to be replicated in the direction of panel width. The data set STRES TABLE contains a 
line of N x , Ny, and N^ values for each different plate wall in the panel cross section. In 
this case, the skin is assigned section property 1, and the blade is assigned property 2. 

The purpose of data set JLOC REPE is to set up a list of local y, z coordinate values 
for a set of joints in the repeating element cross section. The seven joints used to define the 
finite element representation of the repeating element are shown in figure A2. In generating 
this table for example 1, heavy use is made of the EAL "register" capability to allow for 
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simple changes to the dimensions of the panel cross section. The data set EDEF REPE 
defines the element connectivities and section properties for the repeating element. The 
joint numbers used to define the element connectivities are the same local numbers defined 
in JLOC REPE. The six elements defined in EDEF REPE are designated by the circled num- 
bers in figure A2. In EDEF REPE there are four integer entries for each defined element. 
The first two entries are local node numbers for the repeating element. These two numbers 
define a line and, together with NMNX, are used to generate a row (in the x-direction) of 
E43 elements. The third integer points to the material property data in data set MATC BTAB 
for this element, and the fourth integer points to the wall section property data in data set 
SA BTAB. 


Runstream PANEL, GEOM 

The runstream PANEL, GEOM (fig. A3) is the driver runstream for constructing the 
panel geometry information. It calls the general runstream PANEL, JLOC (fig. A4), which 
uses primarily the data set JLOC REPE to produce joint location data for all repeating ele- 
ments in the panel. PANEL, GEOM also calls runstream PANEL, ELD (fig. A5), which uses 
primarily the data set EDEF REPE to define all the two-dimensional E43 elements in the 
panel. Boundary condition information is defined in PANEL, GEOM with subprocessor CON 
in processor TAB. The commands under CON 1 define classical simple support boundary 
conditions for the nodes on all four panel edges including the stiffeners. Material property 
and wall section property information can also be defined in PANEL, GEOM, as is done here. 
For this particular example, since it contains composite laminates, the wall section proper- 
ties are defined in a separate runstream SA, BLADE (fig. A6). Thus, the SA section property 
data produced by runstream SA, BLADE supersede that produced in PANEL, GEOM. 

Runstream PANEL, BUCK 

The runstream PANEL, BUCK (fig. A7) actually performs the buckling analysis of the 
panel defined with the previously described runstreams. It calls runstream PANEL, LOADS 
(fig. A8), which uses the data sets EDEF REPE and STRES TABLE to directly apply the pre- 
buckling stress state to every E43 element in the panel. Then, linear and geometric stiffness 
matrices are produced and the eigenvalue problem solved. As shown, PANEL, BUCK is 
written to solve a shifted eigenvalue problem based on the value of register SHIFT. SHIFT, 
used to speed convergence, is merely a guess at the lowest eigenvalue and can be set to zero. 
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*(28, PANEL, BLADE) ABC 
*XQT A US 
$ 

$ RUNSTREAM THAT PREPARES GEOMETRY DESCRIPTION AND PRE-BUCKLING 
$ STRESS RESULTANT TABLES FOR EXAMPLE 1. 

$ 

! NMNX=3 7 
! NREP=6 
! LPAN=30. 

! B=5.0 

! D=1.352 

! SHIFT=.93 

TABLE (NI=3,NJ=2) : STRES TABLE 
J-l : -811.35 0. 0. 

J-2 : -697.67 0. 0. 

TABLE (NI=2,NJ=7) : JLOC REPE 
! B02=B/2. 

! B04=B/4. 

! B034=3.*B04 

! D02=D/2.0 

J=1 ,7 

0 . 0 . $ 1 

" B04" 0. $ 2 

” B02 " 0. $ 3 

"B02" "D02 " $ 4 
" B02 “ "D” $ 5 
"B034" 0. $ 6 

"B" 0. $ 7 

TABLE (NI=4,NJ=6,RM0DE=2, TYPE-0) : EDEF REPE 
J=1 ,6 
12 11 

2 3 11 

3 4 12 

4 5 12 
3 6 11 
6 7 11 
*RETURN 
* ABC 


Figure Al.- Listing of runstream PANEL, BLADE. 
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Figure A2 . - EAL repeating element for example 1 with 
local joint and element numbers. 


*(28, PANEL, GEOM) ABC 
*XQT DCU 
$ 

$ GENERAL RUNSTREAM THAT CREATES THE PANEL GEOMETRY 

$ 

PRINT 1 JLOC REPE 
PRINT 1 EDEF REPE 
PRINT 1 STRESS TABLE 
*CALL( PANEL, JLOC) 

*XQT TAB 

MATC : 1 10.5+6 .32 .101 
SA 

1 .084 

2 .058 

3 .058 

4 1.0 
CON 1 

! MM1 = NUMN - NMNX + 1 
ZERO 2,3,4 : 1, "MM1" "NMNX" 

ZERO 2,3,4 : "NMNX" "NUMN" "NMNX" 

ZERO 1,3,5 : "MM1" "NUMN" $ SIMPLE SUPPORT 
ZERO 1,3,5 : 1,"NMNX" $ SIMPLE SUPPORT 
*CALL( PANEL, ELD) 

*RETURN 
* ABC 


Figure A3-- Listing of runstream PANEL, GEOM. 


APPENDIX 


*(28, PANEL, JLOC) ABC 
*XQTC U1 
$ 

$ A GENERAL RUNSTREAM TO CREATE JOINT LOCATION DATA FOR A PRISMATIC 
$ STRUCTURE COMPOSED OF REPEATING ELEMENTS. 

$ 

$ INPUT DATA SET : 

$ JLOC REPE 1 1 - TABLE DEFINING NODES IN REPEATING ELEMENT CROSS SECTION 

$ 

*REGISTER EXCEPTIONS* NMNX, NREP, LPAN 
* REGISTER STORE (15 HOLD REG 11) 

$ NMNX - NUMBER OF NODES IN LENGTH DIRECTION 
$ NREP - NUMBER OF REPEATING ELEMENTS 
$ LPAN - LENGTH OF THE PRISMATIC STRUCTURE 
$ 

! NNRE* TOC, NJ(1, JLOC, REPE, 1,1) $ NUMBER OF NODES IN REPEATING CROSS SECTION 
! NUMN* NNRE-1 : !NUMN“ NUMN*NMNX*NREP + NMNX 
! NNRE : ! NUMN 
*XQT TAB 
START "NUMN" 

JOINT LOCATIONS 

! YBASE=0.0 : ! DELY=DS, 1 , "NNRE" ,1(1 , JLOC, REPE ,1 , 1) 

! JBASE=0 : ! DELN=NNRE*NMNX-NMNX 

$ 

$ LOOP OVER ALL REPEATING ELEMENTS 

$ 

! CNT1=NREP 
*LABEL RELOOP 
$ 

$ LOOP OVER THE CROSS SECTION NODES IN THE REPEATING ELEMENT 

$ 

! CNT2=NNRE 
! JJ-1 
* LABEL NDLOOP 

! Y=DS,1,"JJ", 1(1, JLOC, REPE, 1,1) : ! Y= Y+YBASE 

! Z=DS , 2 , " J J" ,1(1 , JLOC , REPE ,1,1) 

! JNT* JJ *NMNX - NMNX + 1 + JBASE 

$ 

"JNT" 0. "Y" ”Z” "LPAN" "Y" "Z" "NMNX” 

$ 

! JJ=JJ+1 

* JGZ , -1 ( CNT2 , NDLOOP ) 

$ 

! JBASE=JBASE+DELN 
! YBASE=YBASE+DELY 
* JGZ , -1 ( CNT1 , RELOOP) 

*XQT U1 

* REGISTER EXCEPTIONS* NUMN 
♦REGISTER RETRIEVE (15 HOLD REG 1 1) 

♦RETURN 
* ABC 


Figure A4.- Listing of runstream PANEL 
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*(28, PANEL, ELD) ABC 
*XQTC U1 
$ 

$ A GENERAL RUNSTREAM TO CREATE ELEMENT DEFINITION DATA FOR A PRISMATIC 
$ STRUCTURE COMPOSED OF- REPEATING ELEMENTS. 

$ 

$ INPUT DATA SET : 

$ EDEF REPE 1 1 - TABLE DEFINING ELEMENT CONNECTIVITY IN THE 
$ CROSS SECTION OF THE REPEATING ELEMENT 
♦REGISTER EXCEPT IONS=NREP, NMNX 
♦REGISTER STORE (15 HOLD REG 1 1) 

$ NREP - NUMBER OF REPEATING ELEMENTS 
$ NMNX - NUMBER OF NODES IN THE LENGTH DIRECTION 
$ 

*XQT ELD 
E43 

! NNRE=TOC,NJ(l,JLOC, REPE, 1,1) 

! NERE=TOC,NJ(l, EDEF, REPE, 1,1) $ NUMBER OF ELEMENTS IN REPEATING CROSS SECTION 
! NELX=NMNX-1 

! JBASE=0 : ! DELN=NNRE *NMNX-NMNX 

$ 

$ LOOP OVER ALL REPEATING ELEMENTS IN THE STRUCTURE 

$ 

! CNT1=NREP 
♦LABEL RELOOP 
$ 

$ LOOP OVER ALL ELEMENTS IN THE REPEATING ELEMENT 

$ 

! CNT2=NERE 
! JJ=1 

♦LABEL NELOOP 

! Jl= DS,1,"JJ", 1(1, EDEF, REPE, 1,1) 

! J2= DS, 2, "JJ" ,1(1, EDEF, REPE, 1,1) 

! Nl= J1*NMNX - NMNX + 1 + JBASE 

! N4= J2*NMNX - NMNX + 1 + JBASE 

! N2= Nl+1 

! N3= N4+1 

! NMAT = DS, 3, "JJ", 1(1, EDEF, REPE, 1,1) 

! NSECT = DS, 4, "JJ", 1(1, EDEF, REPE, 1,1) 

$ 

NMAT=” NMAT" : NSECT=” NSECT" 

$ 

"Nl" "N2 " "N3 " "N4" 1 "NELX" 

$ 

! JJ=JJ+1 

* JGZ , -1 ( CNT2 , NELOOP ) 

! JBASE=JBASE+DELN 
♦JGZ , -1 ( CNT1 , RELOOP ) 

♦XQT U1 

♦REGISTER RETRIEVE (15 HOLD REG 1 1) 

♦RETURN 
* ABC 


Figure A5 . - Listing of runstream PANEL, ELD. 


APPENDIX 


*(28, SA, BLADE) ABC 
*XQT AUS 
$ 

$ RUNSTREAM THAT CREATES LAMINATE PROPERTIES FOR EXAMPLE 1. 

$ 

TABLE (NI=9,NJ=3) : OMA DATA 0 0 

J=1 : 19. +6 1.89+6 .38 .93+6 .0055 0. 0. 0. 0. $ 

J=2 : 19. +6 1.89+6 .38 .93+6 .0495 0. 0. 0. 0. $ 

J=3 : 19.+6 1.89+6 .38 .93+6 .0110 0. 0. 0. 0. $ 

TABLE(NI=2 ,NJ=12) LAM OMA 1 0 $ PANEL SKIN 
J=1 : 1. 45. 

J=2 : 1. -45. 

J=3 : 1. -45. 

J=4 : 1. 45. 

J=5 : 1. 0. 

J=6 : 2. 90. 

J=7 : 2. 90. 

J=8 : 1. 0. 

J=9 : 1. 45. 

J=10 : 1. -45. 

J=ll : 1. -45. 

J=12 : 1. 45. 

TABLE(NI=2,NJ=10) : LAM OMA 2 0 $ BLADE STIFFENERS 
J=1 : 1. 45. 

J=2 : 1. -45. 

J=3 : 1. -45. 

J=4 : 1. 45. 

J=5 : 3. 0. 

J=6 : 3. 0. 

J=7 : 1. 45. 

J=8 : 1. -45. 

J=9 : 1. -45. 

J=10 : 1. 45. 

*XQT TAB 

S A( 0MAL=1 , MAXNL=1 2 ) 

*RETURN 
* ABC 


Figure A6.- Listing of runstream SA, BLADE 
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*(28, PANEL, BUCK) ABC 

$ 

$ GENERAL RUNSTREAM TO FORM STIFFNESS MATRICES AND SOLVE THE 
$ BUCKLING EIGENVALUE PROBLEM. 

$ 

*XQT E 
*XQT EKS 

* CALL (PANEL, LOADS) 

*XQT DCU 
TOC 1 
*XQT TAN 
*XQT K 
*XQT KG 
*XQT AUS 

KSHFT=SUM(K," SHIFT" KG) 

*XQT RSI 
RESET K=KSHFT 
*XQT EIG 

RESET PROB=BUCK,NREQ=l ,INIT=4 
RESET SHIFT=" SHIFT" ,K=KSHFT 
*XQT DCU 
TOC 1 

COPY 1,2 BUCK EVAL 
COPY 1,2 BUCK MODE 
TOC 2 
*RETURN 
* ABC 


Figure A7-- Listing of runstream PANEL, BUCK. 
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*(28, PANEL, LOADS) ABC 
*XQTC U1 
$ 

$ A GENERAL RUNSTREAM TO DIRECTLY APPLY A SET OF SPECIFIED 
$ PRE-BUCKLING STRESS RESULTANTS TO A PRISMATIC STRUCTURE 
$ COMPOSED OF REPEATING ELEMENTS. 

$ 

$ INPUT DATA SETS : 

$ EDEF REPE 1 1 - TABLE DEFINING ELEMENT CONNECTIVITY IN THE 
$ CROSS SECTION OF THE REPEATING ELEMENT 

$ STRESS TABLE 1 1 - TABLE DEFINING NX, NY, NXY VALUES FOR EACH 
$ SA TYPE IN THE REPEATING ELEMENT 
*REGISTER EXCEPTIONS=NREP, NMNX, NUMN 
*REGISTER STORE (15 HOLD REG 1 I) 

$ NREP - NUMBER OF REPEATING ELEMENTS 

$ NMNX - NUMBER OF NODES IN THE LENGTH DIRECTION 

$ 

! NERE=TOC,NJ(l, EDEF, REPE, 1,1) 

! NELX=NMNX-1 
*XQT AUS 

SYSVEC : STAT DISP 

1=1, 2, 3, 4, 5, 6 : J=1,"NUMN" : 0. 0. 0. 0. 0. 0. 

ELDATA : DISL E43 1 1 

$ 

$ LOOP OVER ALL REPEATING ELEMENTS IN THE STRUCTURE 

$ 

! IELT=1 
! CNT1=NREP 
* LABEL RELOOP 
$ 

$ LOOP OVER ALL ELEMENTS IN THE REPEATING ELEMENT 

$ 

! CNT2=NERE 
! JJ=1 
* LABEL NELOOP 

! NSECT=DS , 4 , ” J J " , 1 ( 1 , EDEF , REPE ,1,1) 

! NX=DS , 1 , "NSECT" ,1(1, STRES , TABL ,1,1) 

! NY=DS, 2, "NSECT" ,1(1, STRES, TABL, 1,1) 

! NXY=DS , 3 , " NSECT” , 1 ( 1 , STRES , TABL ,1,1) 

! IEL2 = IELT + NELX - 1 
G=1 : E=" IELT" , " IEL2 " 

"NX” 0. "NY" 0. "NXY" 0. 0. 0. 0. 0. 0. 0. 0. 0. 

! JJ=JJ+1 

! IELT = IELT + NELX 
* JGZ , -1 ( CNT2 , NELOOP ) 

* JGZ , -1 ( CNT1 , RELOOP ) 

*XQT DCU 

CHANGE 1 DISL E43 1 1 KGF E43 1 1 
*XQT GSF 
RESET KGF=1 
*XQT U1 

* REGISTER RETRIEVE (15 HOLD REG 1 1) 

*RETURN 
* ABC 


Figure A8.- Listing of runstream PANEL, LOADS 
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